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Abstract 

Using the theory of Zhang and Clougherty [Phys. Rev. Lett. 108, 173202 (2012); 
arXiv:1012.4405j, we provide detailed supporting information concerning the numerical calcula- 
tions of the probability s (E) for a low-energy electron with incident energy E adsorbing to the 
surface of nanoporous silicon. 



Electronic address: dpc@physics.uvm.edu 



I. INTRODUCTION 



The numerical values for the rate of sticking R calculated in Zhang and Clougherty (ZC) 
[H are obtained by using Eq. (33) in ZC, 
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In this supplementary note, we provide numerical values for the variables in this expression 
for R. 

The coupling constant gi, the strength of bath-assisted particle transitions between the 
continuum and the bound state, is given by Eq. (36) in ZC 
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where Vq{x) is the static surface potential with ground state energy denoted by — E b . The 
coupling constant g 3 is similarly calculated 
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The coupling strength a was introduced in Eq. (5) in ZC 



J(i0) = 22 -ft2° 2 S ( U ~ U q) = aUJ 
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where the frequency of the excitations is cut off from above by u c . The exact form of a (u q ) 
for the case of coupling to Rayleigh phonons is given later in this note (see Eq. (|25|) ). 

Five fundamental variables g\, g^, a, E b and u c are needed in order to calculate the 
numerical sticking rate given by Eq. ([T]). Once the numerical values of these fundamental 
variables are obtained, the other variables are calculated with the following formulae, which 
were also given previously in ZC. 

The constant p is calculated by 
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The variable z is given by Eq. (28) in ZC 
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where K is given by Eq. (29) in ZC 



The probability of sticking to the surface is the sticking rate per surface area per unit 
incoming particle flux, 
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8= y^- R (8) 

In the following sections, we will take the example of an electron sticking to the surface of 
porous silicon (pSi) (porosity P = 92.9%) by emission of a Rayleigh phonon to show how to 
calculate the five fundamental variables for the specific case of an electron sticking to porous 
silicon. The sticking probabilities displayed in Fig. 1 in ZC are calculated using Eq. (8). 



II. FUNDAMENTAL VARIABLES: g u g 3 & E b 



We adopt the commonly used mode 
x c = 0.05 A, plus a repulsive hard wall 

VJx) = 



of an attractive image potential with a cut off at 
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where 
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and q e is the charge of an electron, eo is the vacuum permittivity and k is the relative 
dielectric constant. The average relative dielectric constant for pSi is found from the 
following three equations J4], 
Bruggeman approximation 



Maxwell and Garnett approximation 
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Looyenga approximation 

i4n = (l-P)Kii + P"& (I 3 ) 
where the relative dielectric constant of bulk crystalline silicon (cSi) and air is respectively 
K cSi = 11.68 and n a = 1. We find for pSi (P = 92.9%), k w 1.2 
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There are exact analytical forms for the continuum and bound state wave functions in 
Vq(x). The relevant matrix elements are given by 



9i 
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where 
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Here, ag is Bohr radius and x' c = x c /clb is the lower limit of the integration, c is given by 
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where J\{x) and Nx(x) are the Bessel and Neumann functions of order one, respectively. 
W Kb is the Whittaker function, whose expansion in terms of Laguerre polynomials 

L n (x) for positive non- integer K b is used in the numerical integrations, 

K b {K h - l)e~2 x L n (x) 
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K b is the bound state quantum number and is found from the boundary condition at x = x c 
on the surface 

w*i(^) = » ( 2 °) 

The bound state energies are given by 
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where h is Planck's constant, c is the speed of light and R y is the Rydberg constant. The 
ground state energy in Vq(x) is — E b = —7.76 meV. k is a dimensionless number calculated 



from the continuum incident energy E, 



E = hcRyklAl (22) 
N Kb is the reciprocal normalization constant of the bound state wave function. Thus, 

<=f K^)^' (23) 

We obtain numerical values of the coupling constants gi and by numerical integration of 
Eqs. ( TT4l) and ( IT5l) . We obtain g^ m 1.3 meVA -1 for an electron bound to the surface of pSi 
(P = 92.9%). gi varies with incident energy E and has a range of 57.7 /ieV A -1 < g\< 1.826 
meVA -1 for incident energies in Fig. 1 in ZC. 



III. FUNDAMENTAL VARIABLES: uj c & a 

We consider particle-surface coupling through Rayleigh phonons. The cutoff frequency 
of Rayleigh phonons is approximated by its dispersion relation jsj 

coc = ^CtQc (24) 

where £ is the ratio of speed of Rayleigh waves to transverse waves, q is the transverse 
speed of sound and Q c is the maximum surface wave vector. We take for pSi (P = 92.9%) 

n n 

£ ~ 0.88 and q/q ~ 0.694, which follows from cSi data [6j. (Ref.[7J shows that £ and 
Ct/ci, where q is the longitudinal speed of sound, are largely independent of porosity P.) 
c\ measured from high pSi with P pa 80% is approximately 1680 ms" 1 Assuming q 
and Ct do not vary much for high pSi with P > 80%, we obtain q 1166 ms -1 for pSi 
(P = 92.9%). We take Q c = 27r/a, where a = 5.43 A (the lattice constant of cSi) |8(. Hence 
from equation ( |24l) . we obtain co c ~ 1.19 x 10 13 s _1 . 



The exact form for particle-surface phonon coupling is given in Ref. [5j. For Rayleigh 

waves, 

where Q is the surface wave vector of Rayleigh waves and S is the surface area of the target. 
The shear modulus G and the target's density p f° r pSi (P = 92.9%) are approximated by 
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using Ref. 

Po = Pcsi (1 - P) ~ 166 kg m- 3 (26) 
G = 0.482p c 2 « 230MPa (27) 

We use a mass density of cSi of p c si = 2330 kg m~ 3 jg]. 
The constant -F(cr) is given by |5j 

w ~ 8? lvr=F"7™ ?+ 8 J (28) 



where a is Poisson's ratio, and r is given by 
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For Rayleigh waves, £ is determined by 



(29) 



£ 6 - 8£ 4 + 8£ 2 (3 - 2r 2 ) - 16(1 - r 2 ) = (30) 

Thus for £ = 0.88, we obtain a « 0.03 and F (or) w 0.24. 

Substituting Eq. ([25]) into Eq. (TJJ and taking the continuum limit of the sum over all 
the modes, we obtain for coupling through Rayleigh waves 

where the following vibrational density of states for Rayleigh waves is used 



fl M = S^ (32) 



Using the numbers of g>3, £, q, G, po and -F(cr) obtained above, we find a ~ 0.008 for an 
electron coupled to the surface of pSi through Rayleigh waves. The numerical values for the 
sticking probabilities s(E) contained in Fig. 1 in ZC are now calculated using Eqs. (1) and 
(8). 
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